Superfluidity in a Three-flavor Fermi Gas with SU(3) Symmetry 
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We investigate tlie superfluidity and the associated Nambu-Goldstone modes in a tiiree-flavor 
atomic Fermi gas with SU(3) global symmetry. The s-wave pairing occurs in flavor anti-triplet 
channel due to the Pauli principle, and the superfluid state contains both gapped and gapless 
fermionic excitations. Corresponding to the spontaneous breaking of the SU(3) symmetry to a 
SU(2) symmetry with five broken generators, there are only three Nambu-Goldstone modes, one 
is with linear dispersion law and two are with quadratic dispersion law. The other two expected 
Nambu-Goldstone modes become massive with a mass gap of the order of the fermion energy gap in 
a wide coupling range. The abnormal number of Nambu-Goldstone modes, the quadratic dispersion 
law and the mass gap have significant effect on the low temperature thermodynamics of the matter. 
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I. INTRODUCTION 



The superfluidity in strongly interacting atomic Fermi 
gas and the associated BCS-Bose Einstein condensa- 
tion(BEC) crossover phenomena 1, 2., 3\ have been ob- 
served in experiments|J, IM S B via the method of 
Feshbach resonance. The experimental study of super- 
fluidity in atomic Fermi gas may be important for us 
to understand the solid-state phenomena such as high- 
temperature superconductivity, and may give some clue 
to search for the ground state of the dense quark matter 
and nuclear matter. In the past years, most theoretical 
and experimental studies concentrated on the two-flavor 
systems such as a ^Li gas vi^ith the two lowest hyper- 
fine states (In this paper, we use the word "flavor" in 
particle physics to denote the internal degrees of free- 
dom of the fermionic atoms) . Compared to electrons in 
solids, atomic systems offer more internal degrees of free- 
dom. For alkali atoms, nuclear spin / and electron spin 
5' are combined in a hyperfine state with total angular 
momentum F. While typical electronic systems are con- 
strained to a SU(2) spin rotational symmetry, the total 
angular momentum F can be larger than 1/2, resulting 
in 2F + 1 hyperfine states differing by their azimuthal 
quantum number mp- Therefore, the atomic Fermi gas 
can provide us a way to study the superffuidity with bro- 
ken symmetry higher than the C/(l) one. In this paper, 
we will focus on a three-flavor system with a SU(3) global 
symmetry. Such a system has been investigated in some 
works 8, 9, 10]. 

It is well-known that, associated with the spontaneous 
breaking of a global symmetry, there should be cor- 
responding Nambu-Goldstone(NG) bosons. Such NG- 
bosons dominate the low temperature thermodynamics of 
the system. According to the Goldstone theorem [ill lia | . 
if an internal continuous symmetry group is sponta- 
neously broken down to a subgroup with N broken gen- 
erators, N NG-bosons appear in Lorentz-invariant sys- 
tems, i.e., the number of NG-bosons is equal to the num- 
ber of broken generators. However, from the Nielsen- 
Chadha(NC) theorem|l2, for systems without Lorentz 



invariance the number of NG-bosons can be less than 
the number of broken generators. Let iVi and iV2 be 
the numbers of gapless excitations which have, respec- 
tively, the dispersion laws c<j ~ |p| and w ^^ |p|~ in the 
limit of long wavelength, the number of broken genera- 
tors satisfles the relation A^ < TVi -|- 27V2. For the equal- 
ity between the number of NG-bosons and the number 
of the broken generators, there is an important criterion: 
If {[Qi,Qj]) — for any two broken generators Qi and 
Qj, i,j — 1, 2, ..., A, the number of NG-bosons is equal 
to the number of the broken generators [1^. 

For the three-flavor Fermi gas with SU(3) symmetry we 
will consider in this paper, the ground state of the system 
contains both gapped and gapless fermionic excitations. 
When the SU(3) symmetry is spontaneously broken to 
a SU(2) subgroup with five broken generators, we will 
show with an explicit calculation that there are only three 
NG-modes. Among them, one has linear dispersion law 
and the other two have quadratic dispersion law. The 
reason for the abnormal number of NG-modes and the 
appearance of quadratic dispersion law is found to be 
the fact that, the condition {[Qi, Qj]) = is not satisfied 
due to the density imbalance between the gapped and 
gapless fermions. 

The abnormal number of NG-modes and the non-linear 
dispersion law have been widely discussed in relativistic 
field theory at finite densitvp^ llSl lla| . They were also 
found in the study of two flavor color superconductivity 
in the Nambu-Jona-Lasinio model[l3] where the condi- 
tion {[Qi, Qj]) = is not satisfled due to the lack of color 
neutrality. However, the abnormal number of NG-bosons 
can not be realized in superfluid quark matter and has 
no observable effect, since the color neutrality should be 
imposed via some mechanism such as gluon condensation 
and the NG-bosons should be eaten up by the gluons via 
the Higgs mechanism. In atomic Fermi gas, there is no 
constraint like the color neutrality, and the NG-modes 
are physical degrees of freedom which dominate the low 
temperature thermodynamics of the system. The theo- 
retic prediction of the NG-modes may be tested in future 
experiments via the measurement of the thermodynamic 
quantities. In addition, the mass gap of the two mas- 



sive collective modes found in U3| is very small compared 
with the quark energy gap, while the corresponding mass 
gap in the three-flavor Fermi gas is of the order of the 
fcrmion energy gap, which makes remarkable effect on 
the low temperature thermodynamics. 

The paper is organized as follows. In Sectional we 
set up the model for the three-flavor Fermi gas with 
SU(3) global symmetry. In Section lTTll we investigate the 
ground state of the system and the symmetry breaking. 
In Section nVl we investigate the pair fluctuation around 
the superfluid ground state, identify the NG-modes and 
find their dispersion laws. In Section we discuss 
whether we can recover the five NG-modes. We summa- 
rize in Section IVII The natural unit oi c = h = ks = i 
is used through the paper. 



where Tq — I3 is the generator of the U{1) group and 
Ta — ^{a — 1,2, ...,8) are the generators of the SU(3) 
group. Due to the above symmetry, the system possesses 
nine conserved charges or generators Qa{a = 0, 1, 2, ..., 8) 
given by 



Qa^ Jd^^^J^Ta^. 



(5) 



Like the two-flavor or f/(l) system, for attractive cou- 
pling g we can perform an exact Stratonovich-Hubbard 
transformation to introduce the pair fields 

$/ -- ■^ii^aieapii'fs), $} ^ Ij^iVj^apli)}) (6) 

for 7 = 1,2, 3. With the Nambu-Gorkov fields defined as 



II. THE MODEL 

The physical system we are interested in in this paper 
is an idea system composed of three flavors of fermions 
with attractive interaction. Such a system can be realized 
in cold atomic Fermi gas such as a ^ Li or ^°K gas where 
the three flavors come from three degenerate hyperfine 
states [13 • Generally, the system can be modelled by the 
Lagrangian density 

c^^Uidt + ^ + iA^ + c^nu (1) 

\ 2m / 

where ■0 = (V'l: '02, "03)"^ ^-nd V'^ = (V"!: V'2: V'3) are the 
three-component fermion fields, /^ is their chemical po- 
tential, and m is their mass. We have assumed that the 
chemical potentials of three flavors are the same due to 
the chemical equilibrium. 

It is generally believed that the s-wave channel is the 
dominant pairing channel. According to the Pauli prin- 
ciple, the total wave function of the Cooper pair must be 
anti-symmetric. Using the decomposition 3 (8) 3 = 3 © 6 
for the SU(3) group, the s-wave pairing must be associ- 
ated with the anti-triplet channel in flavor space. In this 
paper we consider only the s wave pairing channel, the 
interaction can be modelled by 



9 



l^int = 7(0a«ea/3/V'^)(V'a'«eQ'/3'/V'/3'): 



(2) 



where g is the bare coupling related to the s-wave scat- 
tering length and e^fe is the total anti-symmetric tensor. 
Throughout, summation is implicit over repeated flavor 
index. Note that the interaction Lagrangian can also be 
written as 



A-'Z' 



E 

a=2,5,7 



(VU,0*)(^^AaV), 



(3) 



where \a{a — 1, 2, ..., 8) are the Gell-mann matrices. The 
model Lagrangian has the symmetry SU(3)® U(l), i.e., 
it is invariant under the transformation 



*=(^*) , *^-(^t ^n 



(7) 



the partition function Z of the system can be expressed 
as 



Z= \\d^\d%\d^\\\d^iY^\^ ' I (8) 



with the kernel /C[<i>j, $/] defined as 



^^ + 2;;i + M if^api^i 



ieat3I^*I -Or 



(9) 



in the imaginary time (r = it) formalism of finite tem- 
perature field theory with J^ — Jq dr J d'^x, where /3 is 
the inverse of temperature, /3 = l/T. Integrating out the 
fermionic degrees of freedom, wc obtain 

Z= /"[d$^][d$/]e-^=^^[*^^*^l (10) 

with the effective action 

5e//[<i>l,<i>/] = / ^^ - iTrln/Ci*^,*,]. (11) 
J X 9 ^ 

III. THE GROUND STATE 

At some critical temperature Tc the system should un- 
dergo the phase transition from the normal phase with 
the SU(3) symmetry to the superfluid phase where the 
SU(3) symmetry is spontaneously broken. Since we focus 
on the low temperature region where T <C T^,, the mean 
field approximation or saddle point approximation is be- 
lieved to be a good treatment for the ground state. The 
order parameters which characterize the superfluid phase 
or symmetry broken phase are defined as the expectation 
values of the pair fields 



i, 



-iTa.ee 



V-, a = 0,l,2,...,8, 



(4) 



($? 



($1), 7=1,2,3. 



(12) 



Let us consider the homogeneous and isotropic super- 
fluid state where the order parameters are independent 
of the coordinates. The thermodynamic potential 17 in 
mean field approximation can be expressed as 



n 



1 



— Seffm = A},'^I 



(3V 
A*jAj 
9 



(13) 



-y 



2/3 



(27r)3 



■TTliig-\iu;n,p), 



where V is the volume of the system and Q ^ is the 
inverse of the fcrmion propagator in momentum space: 



g ^(iw„,p) 



ieai3iA*j iu, 






(14) 



with LOn the Matsubara frequency for fermions and S,p — 
p^ /{2m) — ji. A straightforward algebra shows that the 
determinate of ^^^ in the Nambu-Gorkov® flavor space 
reads 



detg-' = [{^u;^r~ep][i^u^nr~ep-A^y 



(15) 



which indicates that the thermodynamic potential de- 
pends only on the quantity A^ deflned as 



A^^IAil 



IA2P 



IA3P. 



(16) 



The fermionic excitation spectra can be determined by 
detQ^^ = 0. For positive chemical potential, the super- 
fluid phase contains both gapped and gapless fermionic 
excitations, 



C.i,2(p) = ±./C2 + A2, C.3(p)=ep- 



(17) 



With the quasiparticle dispersions, the thermodynamic 
potential can be evaluated as 



n 



mA 



dh 



1 






1 
Ye. 



(18) 



where we have defined the notation Ep ~ ^/Q + A^ and 

replaced the bare coupling g by the low energy limit of 
the two body T-matrix 



47rac 



d^p 1 

l2^2^ 



(19) 



with Qs the s-wave scattering length and ep = p'^ /{2m). 
Now we discuss the symmetry breaking pattern. Simi- 
lar to the two flavor color superconductivity, the symme- 
try breaking pattern in the current case is 



SU{3) (g> U{1) — > SU{2) (g) U{1). 



(20) 



To see the broken and unbroken symmetry groups explic- 
itly, it is convenient for us to choose 



Ai 



0, A3 = A ^ 



(21) 



without loss of generality. In this case, only flavors 
1 and 2 participate in the Cooper pairing and flavor 
3 remains unpaired, a SU(2) subgroup with generators 
Ti,T2,T^ and a t^(l) subgroup with generator Tq = 
(voTq — Ts)/2 remain unbroken, and the broken genera- 
tors are T4,, T5, Tg, T7 and fg = {To+ViTs)/2. We should 
emphasis that all the physical results do not depend on 
the speciflc choice of symmetry breaking direction due to 
the fact that the Lagrangian is invariant under the SU(3) 
transformation. 

To determine physical quantities in the superfluid 
state, we should solve the gap equation together with the 
number equation. Assuming the total number density n 
is fixed, we can introduce the Fermi momentum pp and 
Fermi energy ep through the definitions n ~ 3 xp^/(67r^) 
and ep = p|^/(2to). At mean field level, the gap equation 
which determines the energy gap A can be derived via 
dfl/dA ~ 0, namely. 



lA 



47ra, 



= A 



_fp_ 

(27r)3 



1 - 2f{Ep 

2 En 



1 

2^ 



(22) 



and the number equation can be derived via n 
—dfl/diJ,, namely. 



_fp_ 

(2^)3 



1 



Er. 



'E, 



fiEp) + fi^p) 



(23) 



where f{x) is the Fermi-Dirac distribution function. For 
the specific choice of symmetry breaking direction H21|l . 
the number densities ni, 712 for the paired flavors and n^ 
for the unpaired flavor can be evaluated as 



ni = ^2 



"-3 



(27r)^ 



d3p 

^2^ 
■fi^p) 



E„ 



^fiEp) 



(24) 



which satisfy n = ni + n2 + n^- Once the Cooper pair- 
ing occurs, the number of the paired fermions becomes 
different from the number of the unpaired fermions, 
rii — n2 ^ ns- This difference can be parameterized 
by the ratio a deflned as 



"-1 + "-2 - 2n3 



(25) 



Solving the coupled set of gap equation and number 
equation, we can obtain the gap A, the chemical po- 
tential /i and the ratio a as functions of the coupling 
{pFas)~^- Before the detailed numerical calculation, we 
give a qualitative estimate at T = 0. In the BCS limit, 
the energy gap A is very small and the chemical poten- 
tial is approximately the Fermi energy, the ratio a will 
be very small. In the EEC limit, the chemical potential 
becomes negative and the unpaired fermions disappear, 
the ratio approaches to the limit a — 1. In Fig^ we 
show the numerical results of the chemical potential /i 
and the ratio a at T = 0. The ratio a is always positive 
which means ni = ^2 > ^3. 



What does the nonzero ratio a mean? To answer this 
question, we calculate the expectation values of the gen- 
erators Qa{a — 1,2,. ..,8). They can be calculated via 
the formula 



(Qa) 



-y 



W) 



rTr 



Tj,®Tag{iu}n,V) 



(26) 



where t^ is the third Pauli matrix in the Nambu-Gorkov 
space. The explicit form of the fcrmion propagator in the 
Nambu-Gorkov(X)flavor space can be evaluated as 



t/(iw„,p) 







V 



with the nonzero matrix elements defined as 

-.A _ '*'^n T Cp -n 1 









-iT \ 


g^ 


iT 





go 








-iT 


gf 











af 








GlJ 



(27) 






A 






± 



iuJn ± ^p ' 



(tLo^r - El ■ 

After a straightforward matrix algebra, we find 

(ga)-O, a -1,2,..., 7, 
(Qs) _ ni + n2 - 2n3 _ an 
~^ ~ 72, ~ V2,' 



(28) 



(29) 



According to the commutation relation of SU(3) group, 
we have 

([Q4, O5]) - ([Qe, O7]) = *V3(Q8) - ianV ^ 0. (30) 

Therefore, the nonzero ratio a in the superfluid phase 
means that, the condition {[Qi^Qj]) = which is suffi- 
cient for the equality between the number of NG-bosons 
and the number of the broken generators is not satisfied 
in such a system. However, we now can not conclude that 
the number of NG-bosons is not equal to the number of 
broken generators. 



IV. THE NAMBU-GOLDSTONE MODES 

We investigate now the pair fiuctuations around the 
superfiuid state and examine whether there are five NG- 
modes. If there exist five NG-modes, they must be the 
collective modes associated with the pair fluctuations. 
After the field shift $3 -^ $3 + A, the effective action 
reads 



Seff^ /^-iTrln[g-i+S(<i>l,$,)] 
Jx 9 ^ 



(31) 



0.5 




\^/e. 



-0.5 



-1 -0.5 0.5 1 

(Ppa^)-^ 

FIG. 1: The chemical potential /i scaled by the Fermi energy 
€f and the ratio a as functions of {pFas)~^ at T = 0. 



with the matrix S defined as 










f 








Z$3 


-i$2 \ 










-Z$3 





i^i 


E = 








i$3 -i$2 


Z<i>2 













-z$; 


i<i>l 













V in 


-i<i>l 








/ 



(32) 



To identity the existence of NG-modes and study their 
dispersion laws, we need to evaluate the effective action 
only to the quadratic terms of the pair fields. Using 
the derivative expansion, the effective action up to the 
quadratic terms can be expressed as 



S, 



eff 



$*$ 



jV/ 



1 



9 



Tr[gS($^$,)gS($^$,)]. (33) 



To do calculations in momentum space, we define the 
Fourier transformations 






1 



W 
1 






-^^'n,r+^q■x 



$7(9), 



— 'ii^n,r+«q-x 






^Ka) 



(34) 



for the pair fields, where the four momentum in the 
imaginary time formalism is defined as q = (w„,q) 
with f„ the Matsubara frequency for bosons. We have 
$*(— g) = ($j(g))* due to the above definitions. After 
a straightforward matrix algebra, the effective action of 
the pair fields can be decomposed into three parts: 

^e//[$l,$/] =5i[$t,$i]+52[$2,$2]+^3[*3,$3]. (35) 

The pair fields for / — 1, 2, 3 do not mix with each other, 
which makes the calculation of NG-modes quite easy. 
Each part of the effective action Si takes the form 

1 



[ll]\q)<P}{-q)<Pl{q) + Ilf {q)<P*j{q)<Pi{-q) 
+ Tl]HqWi{-qWi{q)+^f{q)'^i{qni{-q)], (36) 



q 



where the functions ny(g) for / = 1,2 take the same 
form due to the residue SU(2) symmetry, 

uYiq) = nf(-,) = i 

9 



n|^(g) = nf(g) = o, 

and the functions H'^ (q) are given by 



(37) 



nl\q) = nf{-q)^^ + ^Y.[stip)S^ip + <i)], 

nl-'iq) ^ Ufiq) - ^ E l^iP)^iP + 1)] ■ 

p 



(38) 



The functions Hj for I = 1,2,3 are evaluated in Ap- 
pendixes A and B. To identify the existence of the NG- 
modes and determine their dispersion laws, it is conve- 
nient to use the real and imaginary parts of the complex 
pair fields defined as 

^i{x) = ((pi{x) + i<l)i{x)) /V2, 

^*i{x) = i^i{x)~i(j)i{x))/V2 (39) 

for / = 1,2,3. The dispersion laws are determined by 
the zeros of the determinate of the matrix 11^ , 

IlYiq)Uf{q)^Uf{q)Uf{q)=0. (40) 

A. The / = 1, 2 or T4, Ts, Te, T7 Sector 

For / = 1,2, after the analytical continuation ivn —^ 
qo + is, the function 11^^ can be expressed as 



H/ (go, q) = qoH{qo, q) + J(go, q), 



(41) 



where the functions H and J are even function of q and 
hence depend only on q^, see Appendix A. Firstly, let us 
examine whether there are four gapless NG-modes corre- 
sponding to the broken generators T4, Ts, Tg, T7. To this 
end, we take q^ = to calculate the mass gaps of the 
collective modes. The mass gaps of the collective modes 
are given by the roots of the equation 



q^H{qo,0)H{^qo,0) = 0, 



(42) 



where we have used the fact J{qo,0) = 0. Obviously, 
^Q = is a root which gives two gapless NG-modes. To 
examine whether there exist other two gapless modes, we 
need to check whether the equation 



H{0,O) = 



(43) 



is satisfied in the superfluid phase. It is easy to find the 
interesting relation between H{0,O) and (Qs), 



iJ(0,O) 



ni + n2 — 2/13 
A2 



AW 



(44) 



Since a can not be zero once BCS pairing occurs, as we 
have shown in the last section, we conclude that, there are 
only two gapless NG-modes, and the other two expected 
NG-modes become massive. 

We now calculate the dispersion law of the gapless NG- 
modes and the mass gap of the massive modes. In the low 
energy limit go — > 0, q — > 0, we can expand the functions 
H and J as Taylor series of (go, q) at the point (0, 0) and 
keep only the leading terms. To find the dispersion law 
of gapless NG-modes, we take the expansion 



qoH{0,0) 



q^ d^J{qo,q) 



9q2 



= 0. 



(0,0) 



Due to the relation obtained in Appendix A, 



d'^Jiqo,ci) 



d(f 



--HiO,0), 



(45) 



(46) 



(0,0) 



the NG-modes have a quadratic dispersion law near q^ = 
0, 



qo 



q 

2m 



(47) 



It is very interesting that here the quantity m is just the 
fermion mass. For the massive modes, we try to find the 
zero of the function H{qQ, q). In the low energy limit, we 
take the expansion 



H{0,O) + qo 



dH{qo,q) 



dqo 



q^ d'^H{qQ,q) 



(0,0) 



9q2 



0, 



(0,0) 



which leads to the dispersion law 

qo = mi + , 

2m2 

where the mass gap toi is given by 
idHiqo,q) 



mi = -H{0,0) 



\ dqo 

and the quantity TO2 is defined as 



(48) 



(49) 



(50) 



(0,0)^ 



7712 



dH{qo,q) 



d'H{qo,q) 



(0,0) 



9q^ 



(51) 



(0,0) 



In the above calculations, we have employed the trick 
of Taylor expansion which is valid in the low energy limit. 
In Fig|21we showed the ratio between the mass gap nii 
of the massive collective modes and the energy gap A 
of the fermionic excitations as a function of the coupling 
{pFas)~^ a.t T — 0. In the weak coupling limit, the ratio 
approaches to zero, which reflects the fact that the effect 
of nonzero a and hence nonzero -ff (0, 0) can be neglected. 
In this case, we can approximately find five NG-modes 



with linear dispersion law in the energy and momentum 
region 



ii < go < A, mi < VF\q\ < A, 



(52) 



where vp = pp /'Tn is the Fermi velocity, and our result 
is consistent with the numerical calculation in lOj where 
the authors worked in the weak coupling limit and found 
five NG modes. However, in a wide range of the coupling 
such as —1 < (j)Fas)~^ < 1 shown in Fig|2 the mass gap 
nil is of the order of A and hence the order of T^, which 
means that the effect of nonzero a can not be neglected 
and the mass gap mi becomes important for the low 
temperature thermodynamics. 

We conclude that the abnormal number of NG-modes 
and the mass gap m-i make sense in a wide range of the 
coupling and have significant effect on the low temper- 
ature thermodynamics. This situation is quite different 
from the case of two flavor color superconductivity in the 
Nambu-Jona-Lasinio model where the mass gap of the 
massive modes is very small compared with the energy 
gap of the quarks |17j. 



> 3 




FIG. 2: The ratio between the mass gap mi of the massive col- 
lective modes and the energy gap A as a function of (pFtSs)"^ 
at T = 0. 



at any temperature below Tc- Hence we have proven that 
there must be a gapless NG-mode corresponding to the 
broken generator Tg. 

Next we determine the dispersion law of this NG-mode. 
Similarly, we expand the functions Il'^{qo,(l) as Taylor 
series of (go, q) at the point (0, 0) and keep only the lead- 
ing terms. From the explicit form of the functions Hg-* 
evaluated in Appendix B, the Taylor expansions to the 
lowest order take the following form 

nl\qo,(l) = A + qoB + €1^0/2, 
Ilfiqo,ci)=A-qaB + q'D/2, 
nf (go, q) = -A + q^„C/2 + q^E/2, 
^l\qo, q) = -A + qiC/2 + q^E/2, (55) 

where the coefficients A,B,C,D,E are not explicitly 
shown. The dispersion law of the NG-mode is determined 
by the equation 

{A + qoB + ci^D/2) {A - qoB + q^D/2) 
- {-A + qlC/2 + q^E/2f ^0. (56) 

Keeping the lowest order in qo and q, we obtain a linear 
dispersion law 



qo = Ws|q|- 
The velocity of the NG-mode is given by 



A{D + E) 
B^~AC' 



(57) 



(58) 



Since the functions H'^ take the same form as the ones 
in the two-flavor system with broken C/(l) symmetry, the 
behavior of this NG-mode will be the same as the one in 
the two-flavor svstem [2lll2^ . 



CAN WE RECOVER FIVE NG MODES? 



B. The I = 3 or Ta Sector 

For / = 3, the situation is quite conventional. One can 
check that the functions take the same form as the ones 
in the two-flavor or f/(l) system. Let us first identify 
that there exists a gapless NG-mode corresponding to 
the broken generator Tg. To complete this task we need 
only the explicit form of the functions IC^ at go = q = 
and check the relation 

Hii(0, 0)nf (0, 0) - Hi2(0, 0)n2i(0, 0) = O. (53) 

Using the explicit form of the functions Hg-* and the gap 
equation for A, we obtain the relation 

Hii(0, 0) = nf (0, 0) = -Hi2(0, 0) = -H21(0, 0) (54) 



We have shown that in a three-flavor Fermi gas with 
SU(3) gauge symmetry, there are only three gapless NG- 
modes in the superfluid state. A natural question we may 
ask is the possibility to recover the five NG-modes. In 
this section we will argue that there is no way to obtain 
five NG-modes in such a system. 

Firstly, one may criticize that the abnormal number 
of NG-modes may be due to the specific choice of the 
symmetry breaking direction Ai = A2 = 0, A3 = A. If 
we take the following symmetry breaking direction 



Ai = A2 = A3 = A/V3, 



(59) 



we have automatically rii = n2 = n^ and hence (Qs) — 0, 
and we may expect five NG-modes with this choice. How- 
ever, as we have emphasized, the physical quantities such 
as the ground state and the number and dispersion laws 
of the NG-modes do not depend on the specific choice of 



the symmetry breaking direction. In the symmetric case 
with Ai = A2 = A3, the broken and unbroken genera- 
tors wiU be changed, and correspondingly (Qi), {Q4) and 
(Qe) are nonzero. If one works with this choice, he will 
certainly get the same number and dispersion laws of the 
NG-modes as we have obtained. 

Secondly, we may relax the constraint of equal chem- 
ical potentials for the three flavors. For instance, with 
the choice Ai = A2 = and A3 = A, we can set 
fJ-i — tJ-2 — ^J' for flavors 1 and 2 and fi^ = 11 + fio for 
flavor 3. By requiring ni = n2 = n^ we can guarantee 
that all {Qa) are zero and the condition {[Qi,Qj]) — 
is satisfied for any i and j. However, in this case, the 
chemical potentials must be adjusted, a proper value of 
/io 7^ is needed, and the SU(3) symmetry of the model 
Lagrangian is explicitly broken with the broken gener- 
ators T4, Ts, Tgj TV- An explicit calculation like the one 
in [1^ shows that the four NG-modes corresponding to 
the I = 1,2 sector obtain a mass gap hq. In the BCS 
region where the quantity a in the equal chemical poten- 
tial system is small, the corresponding ^0 is also small. 
Such a phenomenon can be regarded as the spontaneous 
breaking of the approximate symmetry, and the corre- 
sponding collective modes with a small mass gap /zq can 
be considered as pesudo NG-modes like the case of color 
superconductivity in the Nambu-Jona-Lasinio model|l9l|. 
However, since the nonzero /^o explicitly breaks the SU(3) 
symmetry, a specific choice of the symmetry breaking di- 
rection may be dangerous i2u\. In fact, if we calculate 
the susceptibilities 



d^n 



dAi 



d^n 



Ai=A2=0 



dAl 



Ai=A2=0 



(60) 



we find they are negative for /io 7^ 0. This indicates that 
if the numbers of the three flavor are fixed, the specific 
choice of the pairing pattern is forbidden. In this case, 
we should solve the chemical potentials ^1, /j,2, /J.3 and the 



condensates Ai, A2, A3 from a large set of equations. 

In conclusion, we can never find five NG-modes in the 
SU(3) system. Only in the BGS limit, we can obtain 
five approximate NG-modes with linear dispersions via 
neglecting the effect of nonzero a. 



VI. SUMMARY 

We have investigated the superfluidity and the associ- 
ated NG-modes in an atomic Fermi gas with three de- 
generate hyperfine states. In our model, the pairing oc- 
curs in the s-wave and flavor anti-triplet channel, and 
the chemical potentials are constrained to be equal due 
to chemical equilibrium. In the superfluid state, there 
are both gapped and gapless fermionic excitations, i.e., 
paired and unpaired fermions. Only in the BEG region 
where the chemical potential becomes negative, the un- 
paired fermions disappear at zero temperature. Once 
the pairs are condensed, the SU(3) symmetry is spon- 
taneously broken down to a SU(2) subgroup with flve 
broken generators. We showed that there are only three 
NG-modes, the one corresponding to the diagonal gen- 
erator is conventional and has linear dispersion law, and 
the other two have quadratic dispersion law. The addi- 
tional two expected NG-modes obtain a mass gap. While 
the mass gap is very small compared with the energy gap 
A of the fermions in the BCS limit, it is of the order of 
A in a wide range of the coupling (ppag)^^ and can be 
reached in experiments of atomic Fermi gas. As a con- 
sequence, the abnormal number of the NG-modes, the 
quadratic dispersion law and the mass gap have signif- 
icant effect on the low temperature thermodynamics of 
the three-flavor Fermi gas. 
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APPENDIX A: THE FUNCTIONS Hy (q) FOR / = 1,2 

In this Appendix we evaluate the functions Hf (q) for 1 — 1,2. From the relations Hj^(g) = H}^(— g) and H}^(g) = 
Il'j^{q) = 0, we need to evaluate H}^ only. After the Matsubara frequency summation and the analytical continuation 
ii^n ^* 9o + *£ we obtain 



Ki<l) = -+l 



72^ 



1 



fiEp) ~ /(^p-,)^2 



QO ^ Ep — £,p-q 



go + Ep- ^. 



P-9I ..2 
p 

S.P-? 



l-f{Ep)-fi^p+,) ^^, 

qO — Ep — ^p_|_q 



f{Ep)-f{^p+,)^ 



qo 



^p ^^^^ ^p v^v.iiiiv.v^ ^^ yj,p 

phase with A 7^ 0, using the gap equation for A, we can express it as 



£,p+q 

(AI) 



with the coherent coefficients ul and w^ deflned as ui — {1 + ^p/Ep) /2 and vi — {1 — £,p/Ep) /2. In the superfluid 



Iiyiq)^qoHiq) + Jiq) 



(A2) 



with the functions H{q) and J{q) defined as 



Hiq) = 

■m = 



(27r)3 2Ep 



qo- Ep - C, 



P Sp-9 



qo + Ep- C. 



p ^p-q 



qo - Ep - S,p^ 



qo + Ep — S^p^ 



rf3p 

(2^ 



Cp-e 



P Sp-g 



2_E„ 



1 - /(i?p) - /(Cp-,) /(i?p) - /(Cp-g) 



qo- Ep- ^, 



P Sp-q 



qa+Ep- ^. 



,p-g 



?P - ep+, f 1 - /(i?p) - fi^p+g) f{Ep) - f{^p+g) 



2E„ 



qo - Ep — ^p_| 



qo + Ep — ^p4 



(A3) 



We now Ust some properties of the functions H{q) and J{q) which are useful to determine the dispersion laws. It 
is easy to observe that the functions H and J arc both even functions of q, 



H{qa,q) ^ H{qo,-q), J(<7o,q) = ^('Zo,-q) 
The function H{0, 0) can be expressed as 



"^^■"'-im. 



l-f{Ep)-,f{ip) , f{Ep)-f{ip) 



Ep +Cp 



^p sp 



By using the identities 



1 



1 



H 



1 



1 



H 



Ep Ep + ^p A^ Ep Ep — ^p A^ 
we derive the relation between -ff(0, 0) and (Qs), 

HiO,0) = ^2 = -^2^- 

The derivative of the function H with respect to go at go = 0, q = can be written as 



aH(go,q) 



10 



(0,0) 



d^p 1 



1 - fiEp) - /(Cp) /(i?p) - /(Cp) 



(^P + Cp) 



{Ep - Cp) 



The function J{q) can be expressed as 

-^ If d'^p p ■ q 

2E„ 



2m m J (27rj 



1 - f{Ep) - /(Cp-,) .f(Ep) - /(Cp-,) 



qo — Ep — ^p_q 



9o + -Ep - ^p-g 



1 - fjEp) - f{^p+,) ^ fjEp) - /(gp+,) 



qo — Ep — S,p+q 

from which wc obtain J{qo, 0) = and 



qO+ Ep - £_p+q 



a"J(go,q) 



dqE 







(0,0) 



for any integer n. For the derivative with respect to q, only the second derivative is nonzero, 

= ^(U,Uj, -—^^ 

(0,0) "^ o^q (0,0) 



dq^ 



0, n^2. 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



(All) 



APPENDIX B: THE FUNCTIONS U^ (q) FOR 1 = 3 

In this Appendix we evaluate the functions 11/' (g) for / = 3. From the relations li^'^{q) — Ilj^{—q) and 113^(5) — 
nP(g), we need to evaluate Ilg^ and Ilg^ only. Completing the Matsubara frequency summation and performing a 



shifting p -^ p — q/2, we obtain 
1 /■ d^p 



nl'iq) = 



"3^(^) -^iw 




9 9 



%^q/2'^i+q/2 '^p-g/2^p+g/2 ^ , , . , .n 

p p „ , p —^ 11 - J{^p-q/2) - J[J^p+q/2) 

90 — ^p-q/2 - i^p+q/2 90 + ^p-q/2 + i^p+q/2 I 



^? -/oW? 






'p-g/2"p+g/2 '^p-q/2'"p+q/2 \ , , . , .n 

——^ '^ ~ W Tp (n^p-q/2) - n^p+q/2)) 

qO + i^p-q/2 — ^p+q/2 QO — i^p-q/2 + ^p+q/2 J 

1 \ 1 - /(gp-g/2) - I{Ep+q/2) 

'^■Ep~q/2Ep+q/2 



1 



90 ^ Ep_q/2 ^ Ep^q/2 % + ^p-9/2 + Ep^g/2 , 



f{Ep-q/2) - f{Ep+q/2) 



1 1 \ f{Ep-q/2) - f{Ep+ 

go + Ep_g/2 — Ep^g/2 go - Ep_q/2 + Ep+q/2 J 2Ep_q/2Ep+q/2 



(Bl) 
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